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Abstract: Structure function of e + e~ —> hadrons cross section proportional to the lon- 
gitudinal part of the hadron tensor is power suppressed with respect to an event shape 
variable in the two-jet region. In the SCET framework, we study the event shape distribu- 
tion for this structure function to NLL level of accuracy. As, a result we obtain the angular 
distribution of hadron jets as a function of the thrust, in the two jet region. We further 
examine effects of non-perturbative hadronization corrections by adopting a shape func- 
tion that reproduce the observed event shape distributions. Impacts of our findings on the 
electroweak measurements via the jet angular forward-backward asymmetry are discussed. 
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1. Introduction 



An intriguing feature of the Standard Model (SM) is the unification of weak and electromag- 
netic interactions. One of the most important predictions of the SM is that the interactions 
of all electroweak gauge bosons are determined by the electromagnetic coupling constant 
a and one additional parameter - the weak mixing parameter sin 2 9\v- 

Study of the Z boson pole at e + e~ colliders provides the most accurate determination 
of the electroweak interaction parameters. The average value of sin 2 6\y found in experi- 
ments at LEP is 0.23153 ± 0.00016 (for a review and bibliography the reader is referred to 
|l[ ) . An experimental method to measure sin 2 6y/ is based on the fact that weak isospins 
for left-handed and right-handed fermions are not the same. This difference of coupling 
constants leads to various angular and polarization asymmetries. The most accurate mea- 
surement of sin 2 9w comes from the forward-backward asymmetry with flavor tagging of 
the final-state quark. The plane orthogonal to the colliding beam lines divides the space 
of all direction into two hemispheres. The electron beam is pointing to the "forward" 
hemisphere. The asymmetry is defined as the relative difference in the numbers of events 
with a reference direction in the "forward" or "backward" hemispheres. Usually the thrust 
axis supplemented with the direction and charge of a prompt lepton from charm or bottom 
meson semi-leptonic decay is adopted as the reference direction. 

Extraction of sin 2 &w with high precision requires taking into account a large number of 
processes accompanying a basic e + e~ — > Z — >• qq process. One of them is QCD interaction 
in the final state. In order to reduce QCD radiation of additional partons, one has to 
implement experimental cuts which effectively select events in the two-jet region 0. The 
event-shape variable thrust (T) can be used to establish such cuts ||. Although, from 
theoretical point of view the thrust distribution is known to unprecedented level of accuracy 
Q, up to now the study of the angular distribution of the thrust axis depending on thrust 
value has not been done beyond the first nontrivial order of perturbation theory ||. It 
was noted in Ref.J|] that an additional contribution to the event shape, which changes 
the angular distribution, is power suppressed in the T — > 1 limit. 

The reason why the thrust axis is so stable in the two-jet region is the following: by 
definition, the thrust axis coincides with the total momentum of final state hadrons in a 
certain hemisphere, thus the multiple branchings, leaving secondary hadrons in the same 
hemisphere, do not change the thrust axis. For the same reason, it is rather hard to estimate 
the influence of multiple hadron radiation on the angular distribution in the strong two-jet 
limit using the present Monte-Carlo event generators. In this paper, we consider in detail 
the angular distribution of the thrust axis depending on thrust value taken in the two-jet 
region. 

The paper is organized as follows. In the next section, we consider possible mechanisms 
for QCD radiations to have an effect on the angular distribution in the two-jet region. In 
SectH using the method of expanding by regions, we find the perturbative correction to 
the structure function, which corresponds to the so-called longitudinal part of the hadronic 
tensor. In sections |I| and |5[ we use the SCET (Soft Collinear Effective Theory) framework 
to study event shape for this structure function, which appears due to a local three-body 
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operator in the effective theory. Comparison with the existing experimental data for the 
angular distribution and the forward-backward asymmetry is left to the last section. 



2. Mechanisms to change the angular distribution 

We restrict our attention to the case where Z-boson is the only intermediate state in 
the processes e + e~ — > hadrons. The thrust (T) dependent cross section for the process 
e + e _ — > hadrons is a contraction of the leptonic tensor and the hadronic one: 



da a 2 TiN c Q 2 

^2 



X _ T dcos0 r dT 2 (Q2_M|) +M|r| 



where Q 2 = (p e - +p e +) , n e = p e -/|p e -|, and ny is the direction of the thrust axis in the 
hemisphere which contains a quark. If one neglects the electroweak radiative corrections 
and the lepton masses, the leptonic tensor reads 

IT (n e ) = (g 2 al + g 2 vl ) 5 f (n e ) - 2ig al g vl a^ (n e ) , (2.2) 

while the hadronic tensor can be parameterized as follows 

= (g 2 vq + g 2 aq ) {F (r) ft (n T ) + 2G (r) gf (n T )} + 2ig vq g aq K (r) a** (n T ) , (2.3) 

where the coupling constants have the following form 

9al = 9ad = -9au = ~ n . , 9vl = 9al (l - 4 sin 2 W ) , 

9vu = ~9ai fl - ^ sin 2 9 W J , g vd = g ai ( 1 - ^ sin 2 %) . (2.4) 



The tensors in Eqs. (|2.2|) and (|2.3|) are defined as follows: 



,r ( u) = + ^ (")<(")+ (2 . 5a) 

(u) = 1 K H " < H [»" ( u ) " »+ H - ( 2 - 5b ) 

a^(u) = ^e^ a V(u)n^ + (u) ; (2.5c) 



where 

2 



u' = l, n=(l,-u), n+ = (l,u). (2.6) 



The expression (2.3) is the most general parametrization of the tensor which satisfies 

(n + n + )^ fl^ = (n + n+) v = 0. (2.7) 

Since we neglect the electron mass and since the hadronic tensor must be contracted with 
the leptonic one, only those structure functions that satisfy the condition (2/7) contribute 
to the cross section even for massive quarks. 
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Let the Z-boson be polarized along the electron beam direction in the e + e~ collision 
rest frame, i.e., it is produced in the the polarized state \l,m) with m = ±1. We then 
obtain the following angular distribution: 



da 



]_ T dTdcos#r 



04 + g 2 aq ) [F (r) (1 + cos 2 9 T ) + 2G (r) sin 2 9 T ] 



m 2g vq g aq K (t) 2 cos 9 T , (2. 



where cos Qt — ii'j 1 ' n e . 

Let us consider a quark-antiquark (qq) pair in the final state in the massless quark 
limit. In the center-of-mass frame, this pair should be produced in the spherical helicity 
state d\ m (n q ) = (1 + \ mcos6 q )/2, where 9 q is the angle between n e and the quark 
momentum p q , and A=±l is the projection of the total spin of the pair on the direction of 
p q . Violation of P-parity by interaction with the virtual boson implies different coupling 
constants for left-handed and right-handed fermions: 



1 1 

9vq ~ 9aq75 = {9vq + 9aq) ^ K 1 ~ 75) + (9vq ~ 9aq) ^ (1 + 75) 

It results in the following angular distribution for the primary gg-pair: 
da 



(2.9) 



d cos l 



(9vq +9aq) \ | +{9vq-9aq) \d_, 



1 



o [04 + 9 2 aq ) (1 + cos 2 9 q ) - m {2g vq9aq ) 2 cos 9 q ] . (2.10) 



Assuming 9 q ~ 9t and comparing Eq.( 2.10| ) with the expression (|2.8| ), we find that the 
following relations hold in the r — > limit: 



F(t) = K(t), 
G(t) = 0. 



(2.11) 
(2.12) 



These relations are the consequence of the free parton approximation and they have the 
same nature as the known Callan-Gross relation || [?]] or the large recoil symmetry relation 
for heavy-to- light form factors 

Additional radiation of high energy partons results in violation of the relations (|2.1l| ) 
and ( p. 12 ). Let us consider radiation of a single gluon with energy E g ~ Q. There are 





(a) (b) 

Figure 1: Three possible directions of the thrust axis in e + e" 
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three different possibilities for the thrust axis to lie along the momenta of the final state 
partons (see Fig.|l|). The topology ((l|c), where the thrust axis aligns the gluon momentum, 
can in principle be excluded from the analysis. To do this, simultaneous tagging of both 
flavored mesons is required. If one tags only one meson then the contribution of the 
topology (|l]c) should also be taken into account. Due to CP-invariance, this topology does 
not contribute to K (r) but gives the main contribution to F (r) and thereby violates the 
relation ( 2.11| ). In contrast, the topologies (Qa, b) contributes mainly to G(r) and hence 
lead to the violation of the relation ( 2.1 2| ) . In the present paper, we mainly consider the 
topologies (|l|a, b) such that \0 q — 9q\ « tt in the r — > limit. We will give a few comments 
about the topology (Jl]c) in Sect.|]. 

In the covariant perturbation theory, the amplitudes are obtained as the sum of the 
two Feynman diagrams, Figs.^(a) and ||(b). For the sake of completeness, we present here 
the leading perturbative results for F (r) — K (r) and G (r) ||: 



tree 



(r) - if, 



tree 



t (I2r 2 + 17r - 45) 



r- 1 



+ | - - 8 In 2 - 4r - 2r 2 J ln(l - 2r) 
+2r(T + 2)lnr + 81n(l-r) [lnr - ln(l - 2r) + 6] + 8 [Li 2 (r) - Li 2 (2r - 1)]} , (2.13) 

V(2-r) 



Gtree (t ) 



a. 



IT 



1-T 



+ 21n(l-r) 



so that the power expansion near r = takes the form: 



Ftree (t) ~ if tree (r) = C F 



7T 



rln- + 0(r 2 ) 



G tree (r) = G(°) (r) + O (r 2 ) = ^C F [r + O (r 2 )] . 



(2.14) 

(2.15) 
(2.16) 



Below we give a qualitative explanation of the physical mechanisms which results in 
violation of the relations (2.11) and ( [2.12 ). Let us consider the the topology (|l]a), where 
the gluon is emitted along the quark momentum direction. In order to remind us of the 
kinematical configuration concerned, we draw the diagrams in Figs. (§a) and (|2|b) in such 
a way that the q, q, g lines follow their momentum directions. 

We find the old-fashioned perturbative theory (oPT) useful to identify two mecha- 
nisms that radiation affects the angular distribution. Using oPT, one can decompose each 




Figure 2: Feynman diagrams for the single gluon radiation process in hadronic Z-boson decays. 
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covariant diagram (Fig. |2|(a) or Fig.^(b)) into the sum of two time-ordered diagrams shown 
in Fig. H and Fig.||[ 

The first mechanism corresponds to the time ordering when additional partons are 
radiated off the primary (/(/-pair after the Z-boson decay, as depicted in Fig.|3[ The in- 
termediate state involved is real but with the different energy Ei, and hance the angular 
distribution ( |2.10| ) is valid for the relative momentum of the primary (/(/-pair. The only 
possibility for the subsequent radiation to change the thrust axis distribution is the follow- 
ing: after the splitting of a parton, one of the "child" -parton is radiated into the "wrong" 
thrust hemisphere, which is opposite to that of its "parent" -parton. For example, the 
oPT-decomposition of the diagram (|2|b) has such intermediate and final states. The main 
effect of such radiation is that the direction of the relative momentum of the real primary 
(/(/-pair can differ from the thrust axis. 

The second mechanism is when virtual states appears before the Z-boson decay, picking 
up the Z-boson from the state d,Q m (ny), as depicted in Fig.|I|. 

Let us estimate the contribution of the first mechanism, assuming that 9 q — 9t misfit is 
of order r 1 / 2 . For the sake of brevity, we assume that the Z-boson is produced in the state 
|1, —1). As one can see in Fig.|2|(b), the primary (/(/-pair is produced in states d±! _ 1 (n g ), 
where n q = P2/IP2 1 , while the thrust axis ny is along the momentum of the antiquark 
pi, i.e. ht = — Pi/|pi I - To derive the distribution with respect to 6t, it is necessary to 
express _ 1 (n q ) through d\ _ 1 (ny), that leads to the admixture of do,-i ( n T) state and 
therefore G{r) =/0 as well as F (r) =/K (r) follow: 

2 i 1 1 2 2 i 1 1 2 

do- ~ (g vq + g aq ) |gLi_i (n q )\ + (g vq - g aq ) ^ (n ? )| 

1 

= Tra^Mdl^M, (2.17) 

n,l=— 1 

with 

T n i = (g vq + g aq f D*\ n Dl u + (g vq - g aq f D^ n D{i, (2-18) 
where D{ is the operator of finite rotations: 

Di m = D{ m (a,e qT ,j). (2.19) 

Here, 6 q T is the angle between the vectors nj- and n q , a is the angle of rotation about 
n^-axis transferring the vector n e x ny into nj- x n q , 7 is the angle of rotation about 




Figure 3: Two particle intermediate state Figure 4: Four particle intermediate state 
in oPT . in oPT . 
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rig-axis transferring the vector x n q into n e x n q . All rotations are denned with respect 
to the "right hand grip rule" . The matrix T n i entering Eq. ( |2.17| ) does not depend on 7 and 
becomes diagonal after averaging over directions transverse to the thrust axis: 



2tt 



d0! rp (2,2 

T n i = {g aa + 9 



x diag 



aq 1 zvq) 
2 , 



1 + cos'' Q qT 2g aq g vq 



+ 



9aq 9vq 



■ COS I 



'qT, 



sin 



1 qT, 



1 + cos 2 9 qT 2g aq g vq 



9aq 9vq 



COS I 



(2.20) 



and the angular distribution of the thrust axis becomes 



da~- 



1 + COS QqT ^ + g 2J ^ + cog 2 + cog (2g aq g vg ) (2 COS 6 T ) 



+ sin 2 9 qT - (g 2 aq + g 2 vq ) sin 2 T . (2.21) 



The coefficients before the distributions 1 + cos 2 9t and 2 cos 9t are now different, but this 
difference is strongly suppressed in the collinear limit 6 q T —> 0: 



1 + cos 2 q T 



2 COS q T 



9 4 



l + -f" + 0(«Jr) 



(2.22) 



1 qT- 



Therefore, the corrections ( [2.21| ) integrated over the region 9 2 T ~ r generate the following 
corrections 1 to the relations (|2.11| ), ( |2.12 ): 



F (t)-K(t)~ t 3 , G~t 2 



(2.23) 



To estimate the effect of virtual states depicted in Fig.|3], let us consider the following 
projectors 

j^(u)^ + (u) h t^+(u)^(u) 



P(+) (u) 



(2.24) 



4 4 

where 4- vectors n and n + are defined in Eq.([2.6|). They satisfy the following projective 
relations: 



Pi 



(±) 



(±)> 



P (±) P W = 0, P {+) + P { _) = 1. 



(2.25) 



To ascertain their geometrical meaning, we consider the following light-like 4-momenta 

q = E(l,v), q' = E(l,n), (2.26) 

and a spinor such that = for a definite helicity A. One can show that the wave 
functions Pr±\ (n) ttq A ' correspond to the two states with the projection of the angular 
momentum ±A on the axis n, but with the same helicity: 



ii, 



W = P (+) (n) u%> + (n) <' = d%l (tt 



(A) 



Aiv) tt_ q / + a x x (0 nv J u q , , 



(2.27) 



1 Here, we do not distinguish between t° and r™ln m r, where the logarithmic part can appear due to 
integration over the gluon energy. 
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where cos 6 nv = n • v. Now let us consider the fermion propagator inside the Feynman 
diagram (^a) . The total momentum of the quark and the gluon is 

=P^+P^ = (p-n)^ + (p-n + )—. (2.28) 

Here and below we omit the arguments of and implying (n^) and (nj-). Using 
the projectors ( |2.24 ), one can split the propagator into two parts 



~ p (-)p^a P{+) + P{+) WTlo p ^ (2 - 29) 

which correspond to retarded and advanced propagations: 

P t P _ p- n 4+ _ 1 4+ / 2 30a) 

( ) p 2 + iO (+) p 2 + iO 2 E - Ex + iO 2 ' 1 ' ' 



P( + ) ^—P( ) = 4^ = ± (2.30b) 

(+) p 2 + iO ( ' p 2 + i®2 E - E 2 - iO 2 



where the following simple relations were used: 

p 2 = (p-n + )(p-n), (2.31) 
p ■ n + = p° - |p| = Q - (|pi| + |p|) = E - E x , 

-p. n = -p°-\p\=Q-(Q + |p| + | P2 | + | P3 |) = E - E 2 . (2.32) 

Here Eq = Q is the energy of the initial state, E\ = |pi| + |p| and E% = Q + |p| + | p>2 1 + |ps| 
are the energies of the intermediate states depicted in Figs.(^) and (|j) respectively. Taking 
into account the decomposition ( [2.27 ), one can see that the term ( |2.30a ) corresponds to the 



quark propagating in the direction while the term ( |2.30b| ) corresponds to the antiquark 
propagating in the direction — ny. Therefore, we can conclude that in the first diagram 
(Fig.^) the Z-boson decays into the quark-antiquark pair with opposite helicities, i.e., from 
the d\.^ _ 1 state. In the second diagram (Fig.f|), the Z-boson disappears being absorbed 
by the intermediate antiquark. Taking into account helicity conservation for massless 
quarks, we conclude that it is possible only from the state c?q According to the energy 
conservation for the final state, i.e., |pi| + |p>2 1 + | p>3 1 = Q, we find that Eq — E 2 = —Q, 
i.e., the propagator ( 2.30bj ) does not depend on the kinematic configuration of the final 



state. Since the phase space for a definite thrust value is of order t, the intermediate state 
shown in Fig. || leads to a contribution 

G(t)~t. (2.33) 

This explains the tree- level result ( p.!6[ )). 

The analysis presented above remains valid in the case of a multiparticle final state 
with primary q and q radiated into the opposite hemispheres. A QCD cascade starting 
from the quark-antiquark state as depicted in Fig. [| cannot change the relations ( [2.11 ) and 



( 2.12 ) significantly in the two jet region, where the corrections are suppressed by r and 



t 2 respectively. In the small r region, the leading corrections to the relation ( 2.1 2| ) come 
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from the process at short distances shown in Fig.|3], which gives rise to the contribution 
of order r ( 2,33| ) . The proper consideration of violation of the relations (2.11) and ( 2.12] ) 
requires a consideration of a new type of jets. As it will be demonstrated below, those 
jets are initiated by two collinear partons produced at short distances, i.e., due to a local 
operator with more than two fields, so that all secondary collinear or soft particles are 
radiated coherently by these partons. 



3. Perturbative corrections 



If one has an integral over a scaleless domain whose integrand depends on a small external 
parameter A< 1, then the method of expanding by regions gives an asymptotic expansion 
with respect to A. The series may contain arbitrary non-integer powers of A as well as 
integer powers of In A. Here, we outline the prescription used below without discussing 
details of the method, for which the reader is referred to the orig inal studies g 0, |l|. 

The method utilizes the fact that the expansion of an integrand with respect to A 
may be invalid in a certain region where the integrand becomes singular. Analyzing the 
integrand singularities, one can establish the so-called power counting rules, according to 
which one can pick up a simplified singular behavior of the integrand by expanding it not 
only in the external parameters but also in integration variables. Using this method, one 
can represent the integral as a sum of integrals, such that each integrand is an expansion 
of the original one with respect to the power counting rules. 

Let us demonstrate how to apply this method to calculate the leading perturbative 
correction to the structure function F (r) in a region where 1-T<t<1: 



F(r) 



{V - 2)N C Q' 



f dp x e(^| P/l -n T |-(l-T)Q) 
J \hex J 

xj>f (nr) (X 



X 



Ju 



a 



F{r) = l + ^F {1) (r) + 



I7T l£) 2 ^ (T) + 



0), (3.1) 
(3.2) 



where = ^„7^V'g an d dp x is the phase space of a final state \X): 



d V Ph 



d Px = (2«r8»\Q--^-Y J Pk ) n , o _, r _ 



hex / hex 



(2k) 



T$(Ph)®{Ph-n + p h -n + ) (3.3) 



and c denotes the sum over spin and color states and V = 4 — 2e is the space-time 
dimensionality. 

First of all, we introduce a small parameter A such that A 2 ~ r. We use the Sudakov 
decomposition to represent all real or virtual particle momenta: 



. n + n 

Ph = {Ph-n) — + (p h ■ n+) - + p h ±. 



(3.4) 
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Region 


Scale 


Power counting Q 1 (p h ■ n,p h ±,p h ■ n + ) 


hard 


Q 2 


(1,1,1) 


r-collinear 


tQ 2 


(1,A,A 2 ) 


/-collinear 


tQ 2 


(A 2 ,A,1) 


soft 


r 2 Q 2 


(A 2 , A 2 , A 2 ) 



Table 1: Power counting rules 



Power counting rules should estimate the components (ph • n,ph±,Ph ■ n+) in comparison 
with A. The regions of integrations and the corresponding power counting rules are pre- 
sented in Table [l]. There are two regions, namely, hard and soft ones, where all momentum 
components are of the same order. In the hard region one should expand an integrand 
with respect to r only and integrate over real or virtual particle momenta in dimensional 
regularization, ignoring any soft or collinear infrared singularities. The expansion in the 
soft region corresponds to an integrand approximation near the soft singularities. There 
are also r- and /-collinear regions, where symmetry among the space directions is strongly 
broken. These regions account for the collinear singularities. 

The contribution of each region is gauge invariant by itself. In Feynman gauge, vir- 
tual corrections to the amplitude (qq\J^\0) contribute in the hard region only, and the 
corresponding contribution to F (r) has the form: 



Fffi* (r) = C F ( 2- 



Q 2 



e A e 3 



(3.5) 



In contrast to the virtual corrections, the real emissions give contributions in the collinear 
and soft regions: 

Fj~ co1 (r) = F[- co1 (r) = C F ( ( 4 + - - vr 2 + 7 + O (e)) + 0(A 2 ), (3.6) 



(i) 



111 \ fx 

^f(r) = C F (^ 



e 2 e 



2^2N-e / 4 ^2 



+ ^T + 0(e) \ +0(X 2 ). 



(3.7) 



The sum of all contributions is a well known result [12]: 

F(r) = l + ^C F f-4m 2 - + 6mi-2+^ + O (a s r) + O (a 2 s ) . (3.8) 

4tt \ T T 3 / 

The singularities with respect to e and the /U 2 -dependence drop out of the sum of all the 
contributions (|3.5|), ( |3.6| ) and (p?^), which tests our use of the method of expanding by 
regions. 

Let us apply the same method to calculate the perturbative result for G (r) function 
in the region 1 — T < r <C 1: 



G(t) 



2tt 



dp 



cr.c 



x 



X 



Jll 



0)G 



\hex 



\Ph ■ n T | - (1 - t)Q 



(3.9) 
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where 

i|| (3-10) 

First, we consider the case when the thrust axis is aligned with the antiquark mo- 
mentum p\ = E\n. The diagrams we would like to consider are depicted in Fig.|2|. The 
momentum p\ is Z-collinear, the momenta P2 and p% are r-collinear in accordance with the 
classification of Table |]. The corresponding amplitudes for the diagrams Figs. (|2]a) and 
@}), which contribute to G(t), have the following form: 



u(P2)Vffv (Pi) = -g s t a u{p 2 )r 



h+fa 



(P2+P3) +*0 V 2 

4 4+\ ~{i>i+h) 



(3.11) 



2 2 / (pi + P3) + * 



7 M v (p 



u {pi) V$)V (pi) = -g s t a u (p 2 ) 
respectively. The sum of the amplitudes ( 3.1 1| ) has the following A-expansion 

fife) + 

(P2+P3) 2 



u(p 2 



v(pi 
2g s t a 



Q 



71 



U (pi 



1 + 



Q 2 



P3 • n + iO 

Using the light-cone gauge, such that the gluon propagator has the form: 



+ 0(A 3 ). (3.12) 



d D xe 



-IP3-X 



<0|T^(x)^(0)|0) 



p\ + iO 



p^n" + p^n^ 
P3 • n + iO 



(3.13) 



and neglecting the power suppressed term (p2 + P3) 2 /Q 2 ~ A 2 , we find that the following 
effective vertex 

_ 2 5 S .0 44+^^44+ 

V eS - V {a) V {b) — - E 



"7" 



(3.14) 



Q 4 ' 4 

gives the leading contribution to G (r) . The r- and /-collinear regions give equal contribu- 
tions, and hence we find: 



G (o) 



167rg 2 C F 

Q A 
a s Cp 



-/ \44+,±, \44+ 1 s 



e(rQ 2 -pi) dp 3 



(3.15) 



7T 



r + 0(r 2 ), 



where p| = (p2 + P3) 2 , e (P3) is the gluon polarization 4-vector e (p$) ■ p% = e (ps) ■ n = 0, 
and d/93 is the element of the three particle phase space fl3.3|) . 

The expression ( |3.14 ) remains valid in the light-cone gauge ( 3.13 ) even when the parti- 
cles are off shell (p 2 ~ p\ ~ p| ~ A 2 Q 2 ), i.e., it appears as the internal part of the expansion 
of amplitudes with more than three collinear particles in the final state 2 . 

The perturbative correction to G^ (r) has the form fl3"fl : 



G(t) = G (0) (r) 



1 + (r) 

47T 



(3.16) 



2 See Ref.|tLa| for details. 
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where G^ (r) contains hard, r-collinear, Z-collinear and soft contributions: 

(r) = + G« d (Q 2 ,n 2 ) + G« ol (rQV) + G« col (tQV) + (r W) . 

(3.17) 



where the term — /3o/e appears after a s renormalization in the vertex (|3.14| ). The hard 
contribution consists only of virtual corrections to the amplitudes (Fig.|2|), where the loop 
momentum is hard (Table p]) : 



^(QV) = (> 



AC F 1 



+C A (-16+^ + 16C(3)^) + C f (34 + tt 2 - 32C(3)) + O (e)\ . (3.18) 



As we shall see later, the hard contribution is by construction the matching coefficient 
of the weak current to the effective three-body operator, whose leading matrix element 
contains the vertex ( 3.14| ). 



In contrast to Eq.([0]), Gjp and G^ col contributions are different, because the 
leading amplitude contains two r-collinear particles and only one /-collinear particle. Ad- 
ditional /-collinear particle generates the same correction as Eq.( |3.6[ ) for F(t), except for 
an additional integration over the invariant mass p 2 ^ of the r-collinear particles: 



Q 2 J Q 2 



tQ 2 \ c „ f 4 7 



(3.19) 



The r-collinear contribution includes virtual corrections, so that the loop momentum is 
implied to be r-collinear (Table |l|), as well as real radiation of three r-collinear and one 
/-collinear particles: 



G 



(i) 

r— col 



(tQ 2 ,^ 



V 2 



! ( ' \ - ^ - 22C(3)) + Cf (-f - y + 44C(3)) - | iV f T F + O (e) }. (3.20) 



The soft radiation also has the same form as Eq.(3.7) for F (r), but with an additional 
integration over the invariant mass p 2 ^ of the r-collinear particles: 



2 2 ,m 2 ) 



47T 

a s T jo 
r 2 Q 2 
A* 2 



soft / ^ PR \ d ^R 



Q 2 

8 vr 2 
e 3 



16 



(3.21) 
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The singularities with respect to e and the ^-dependence drop out of the sum of all 
the contributions: 



G« (r) 



1 



1 



4C F ln z - + In - 

r r 



y + 12C(3) 



+ C A 



14 



4 

-T F A^ f 
3 



™T F iV f . 



(3.22) 



It is interesting to compare this correction with that of the leading thrust distribution 
Q3.8| ). The double logarithmic corrections are equal. It has a simple physical interpretation: 
resummation of double logarithms results in a statistical factor responsible for excluding 
some part of radiation which is out of the two-jet region, and hance this factor is insensitive 
to the jet structure. However, the sing le logarithms in Eq.(U) and Eq. (B) are different. 
In the case of F(t), it gives a positive contribution because GCplnr^ 1 > 0, but the 
correction to the G(r) distribution ( 3.22 ) contains the logarithm lnr -1 with a negative 
coefficient: 



4^ 
3 



14 + C A 



23 

y 



2f_ 
3 



-T F N f = 1 



2n 2 



< 0. 



(3.23) 



The nonlogarithmic correction to (r) is positive and larger than Cf (27r 2 /3 — 2) in 
Eq,(|3.g|). We will find in section [5] that these differences affect the event shape distributions 
F(t) and G{t) after resummation of large logarithmic corrections. 



4. Factorization formulae 

In this section we give the heuristic derivation of the factorization formula for G(r), which 
allows us to sum all big logarithmic corrections to the NLL level of accuracy. 

We first note that power counting rules imply a hierarchy of the components of integra- 
tion momenta and hence a hierarchy of the components of fields. Instead of an expansion 
of integrands obtained in a full QFT, it is sometimes possible to introduce an effective the- 
ory such that its Feynman rules reproduce the expanded fundamental amplitudes. Field 
modes corresponding to different regions are associated with different fields of various effec- 
tive theories. Although, in perturbative calculations, the effective theory framework does 
not provide new information, it turns out to be extremely efficient when one needs to es- 
tablish factorization formulae and evolution equations for resummation of large logarithms. 



A lively presentation of these ideas is provided e.g. in Ref.|14], where, starting from the 
minimal set of assumptions about hadronic final states, the authors use the soft collinear 
effective theory (SCET) |L5], [ll| to derive the factorization formula for the distributions of 
a large class of infrared safe observables (angularities). 

In order to construct the SCET Lagrangian for collinear quarks interacting with 
collinear and soft gluons, one has to split the quark field tp using the projectors fl2,24|) : 



i> = in + Vn, in = P(+)*P, Vn 

and split the gluon field into the collinear and soft parts: 

A. = A c ^ n + A s . 



(4.1) 



(4.2) 
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The key assumption is that all the collinear field operators £ n , r\ n and A c%n , A s generate field 
modes with momenta obeying the power counting rules listed in Table |l[ Here collinear 
quanta corresponding to A c ^ n are /-collinear and soft quanta A s have soft momenta. This 
assumption implies that we restrict our consideration to the processes which can be de- 
scribed by the QCD fields that are smooth in certain directions. The power counting rules 
of Table |l| lead to the following hierarchy of the components of fields: 



a, ~ y 



in ~ A, 

n+ ■ A c>n ~ 1, A c>n ± ~ A, n ■ A CiU ~ A 2 

where the Sudakov decomposition of the collinear field A c n is used: 



(4.3) 



n> 



n r . 



(n+ ■ A C)n ) — + (n • A Cjn ) — + A c , n± . 



(4.4) 



Using the estimations ( [4, 3D , one can expand the QCD Lagrangian to all orders in A. This 
expansion truncated to some order in such a way to preserve the invariance under the 
homogeneous gauge transformations [[j/?]] corresponds to the SCET Lagrangian [15, 16]. 
Here, we outline the result of Ref. [O] for the thrust distribution: 

F{t) = H (Q 2 ,fi 2 ) J dp 2 h d P ldk J (plfi 2 ) J {p 2 R , !i 2 ) S T (k, fi 2 ) G (Q 2 t -pl-pl- Qk) . 

(4.5) 

H (Q 2 ,/j?) is the hard function or the square of the usual on-shell QCD Sudakov form 
factor. J (p 2 , /u 2 ) is the jet function: 



1 



(p ■ 71+) iV c 

x — Im 

2vr 



d 4 xe~ tpx ( 



TUnW^w^(o)e;(o) 



(4.6) 



that is, up to an overall factor, the imaginary part of the QCD quark propagator in the 



light-cone gauge {^' n and and Wilson line W n are defined below in Eq.( Ol )). The soft 
factor St (k,/J?) is defined as follows: 



S T (k,v 2 )=^\(x 



x 



1 n + I n 



o 



\ 2 

J S(k-n-px L -n + -px R ). 



W and Y are collinear and soft Wilson lines, respectively: 



W n (x) = P exp 
Y n (x) = P exp 



f 

/ ds n+ ■ A c (x + sn+) 

J —00 



ds n ■ A' (x + sn) 



(4.7) 

(4.8) 
(4.9) 



In the effective theory framework the vertex ( |3.14j ) corresponds to the following oper- 



ator: 



3 = 3R + 3L , 3R = 2g s L4 ± U 3L = 2g s L4 ± , n £ n , (4.10) 
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where £ nj and A± n . are fields from different SCET copies in the light-cone gauges of the 
type of Eq.( |3.13| ) with the light-like vector m. In order to introduce the gauge invariant 
operators, one should replace the fields entering the operators ( |4. 10|) to f and A'j_ in an 
arbitrary gauge, using the following relations: 



£ = YW^', g s A, 



c± 



Y (wUD' c± W 



di)Y\ 



D f cL = d ± - igA[ 



(4.11) 



c± 



where fields without primes are in the corresponding light-cone gauge. Including a soft 
Wilson line Y in the definitions ( |4, 11 ) allows one to decouple soft and collinear degrees 
of freedom in the leading order SCET Lagrangian |l^, [17|. Using the expressions ( 4.11| ) 
yields the following operators 



3 = 2g s i' n+ 4^ n W n+ Y r \ + YnWlCn + 2g s ? n+ W n+ Yi + Y n W n 4 Xt J' n , (4.12) 



where 



A, 



A', 



(4.13) 



The operator (4.12) is in fact the operator 3 derived in Ref. |19|| taken in the limit: n 



n+, Ug 



n and n 



n or n„ 



Integration over hard modes gives the matching 



coefficient C H of the QCD operator J\\ ( pTToj ) onto the SCET operator 3 ( gig ) 

J\\ ^C H (Q 2 ,fi 2 )0 3 . (4.14) 

The important point about the operator (|4.12|) is that it is a local product of the 
r-, /-collinear and soft SCET operators. According to Ref. [14], this feature is the only 
requirement to establish a factorization formula for an angularity distribution. For the 
operator ( [112 ), the thrust distribution takes the form: 

G (r) = 2 H 3 (Q 2 ,^ 2 ) J d P 2 L d P 2 R dk E x (pi, /u 2 ) J (pi, fi 2 ) 

xS T (k, f i 2 )Q(Q 2 r-p 2 -p 2 R -Qk), (4.15) 

where St (k,fi 2 ) is the same soft factor as defined in Q4.7|) , J (pj^/i 2 ) is the jet function 
defined in (p[), and #3 = |Ch| 2 . The new object in the formula ( |4.15| ) is (p^, /u 2 ) : 



x Im 



(p • n) Q 2 N C tt 



d v xe- tpx ( 



n+ 



(*) I (Wl4±,n + C + ) (°) 







(4.16) 



which can be considered as the imaginary part of the quark "transverse" self energy pro- 
jected onto ijt. In contrast to the jet function or the soft factor, whose leading expressions 
are (5-functions, the tree level expression for Sj_ (p 2 ,^ 2 ) is a smooth function: 



,(0) 



Us jy 2 ) C F ( pi 



4ttQ 2 



D/2-2 



2r (D/2) 

r (v - 2) 



(4.17) 
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The most important property of the operator ( 4,12 ) is that it is a local product of 
the collinear fields. The reason is that dominant underlying process shown in Fig. || does 
not depend on the kinematics of the final state since the energy of the intermediate state 
E2 = Q + |p| + I F>2 1 + |k| is equal to 2Q due to energy conservation for the final state. 
It results in locality of the vertex ( p. 14 ) for the three-particle interactions. The essential 



consequence of the locality is that soft gluons are radiated off by just like the two-prong 
QCD antenna ( |4.7| ) for the leading contribution to the structure function F; see Eq.([4.5|). 
An intuitive explanation of this fact is that soft gluons are coherently radiated off all 
collinear emitters and the factorized amplitude of the soft radiation does not depend on 
the fraction of the collinear momentum carried by a particular emitter. For example, using 
a noncovariant gauge, one can show that the amplitude of soft radiation off the r-collinear 
pair depicted in Fig.[| is equivalent to the corresponding amplitude for a single quark 3 : 

t a t t e (*)•«+ _ . fabcf e(k)-n + = ^ e (AQ • n + _ 
k ■ n + k ■ n + k ■ n + 

All those remarkable properties of the structure function G such as the locality of the 
corresponding effective operator and the universal soft factor can be hardly generalized to 
the difference F—K, where the third topology (Fig.|l|c) accounts for the leading contribution 
Q2.15D - In this case the soft radiation undoubtedly differs from that for the topologies 
(Fig.|l|a, b) because one hemisphere contains the single high energy gluon only. Moreover, 
the part of the cross section contributing to F — K is singular in the region where the energy 
of the quark or antiquark tends to zero. This singularity, which is quite similar to that in 
the 7*7 — > 7T° process, leads to an additional lnr in the difference F - K in Eq.( ^15|) . It 



indicates that the effective operator for F — K should be nonlocal in light-cone directions 
and the factorized objects should be different from those used in Eq.( 4.15| ). We will discuss 



F — K contribution to single flavor tag measurements briefly in section |?j. 
5. Resummation of large logarithms 

In order to perform integration in the factorization formula ( }4.15| ) we take the Laplace trans- 
form of each function entering Eq.( |4.15 ) except for the hard coefficient function H3 (Q 2 ,/j, 2 ) , 



G (r) = 2 H 3 (Q 2 ,^ 2 ) ^.\% E ± (sQ 2 ,^ 2 ) j (sQ 2 , /i 2 ) s T (sQ, » 2 ) , (5.1) 

where s = 1/ (^Q 2 e 7E ) and 

j(sQ\» 2 )^ dpV^J(pV), ar(8Q,f)= dke-^ k S T (k,v 2 ), 
Jo Jo 

poo 

[sQ\ V 2 ) = / dp 2 R e~^i S ± {pi, n 2 ) . (5.2) 



Now we use the fact that the expression ( |5.1| ) does not depend on fi 2 . We can exclude all 
collinear logarithms by setting fi 2 = tQ 2 ~ A 2 . In doing so, we can neglect all higher-order 



3 In fact, the same argument is the same as that for explaining the angular ordering in QCD sequential 
branching process (Ref.ffl |1 1|]). 
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corrections in Ej_ (sQ 2 , // 2 ) and j (sQ 2 , // 2 ) since all of them contribute either to the N 2 LL 
level, or to the (pre-exponential) factor, which does not depend on r and can be found 
from the fixed-order result ( p. 22 ), Therefore, we can replace in Eq.( |5.1| ) the jet function 
j (sQ 2 ,fi 2 ) by unity and t ± (sQ 2 ,/x 2 ) by E ± (sQ 2 ,tQ 2 ): 

El (sQ 2 , tQ 2 ) = r dp|e-^ 4 0) (Pi, rQ 2 ) ^ ^ ° F 1 

Jo 



V=4 



2tt 



^Q 2 



(5.3) 



Thus, we obtain the following distribution: 
2a s (tQ 2 ) C f 



G t) 



(5.4) 



7T ' 2ni Jq v 2 Q 2 

Let us now consider the evolution equation for the hard coefficient functions and the 



soft factor | 



dHj(Q 2 ,^ 
din [j? 

ds T {sQ,fi 2 ) 
dln/i 2 



r cusp [a 8 {fi 2 )} ln^J + 7 H ° [a a (/z 2 )] \ H { (Q 2 ,/i 2 ) , 



r cusp [ats (fJ 2 )] In 



s 2 Q 2 



7 S [« s (/i 2 )] Lr(sQV), (5-5) 



where 



r cusp (a s ) = — r (0) + (— j r (1 



47T 

a 



7 l"s 



/a s y 



+ .... 



(5.6) 



In the NLL accuracy, one needs the following expressions: two-loop r cusp , two-loop a s (/i 2 ) , 
one- loop 7*: 

4C F 



MO) 
7(0) 



4C F , 
-6C F , 



(i) = 

7(0) 



9 



[C A (67 - 3ir 2 ) - 20T F N f ] , 



(5.7) 



is found in Ref . , 7^ and 7^ 3 can be found from the expressions ( |3.5D and fl3.1£j ), 
respectively. The initial conditions Hi(Q 2 ,Q 2 ) and st(sQ,sQ) for the equations ( |5jj ) 
contribute to the pre-exponential factor. 

Let us, for a moment, omit all pre-exponential factors and set fi 2 = tQ 2 , thus the 
solutions of the equations (|5.5|) can be represented as follows: 



where 



Hi {Q 2 ,tQ 2 ) = exp {F Hi [L, a s (Q 2 )] } , 
s T (sQ,tQ 2 ) = exp jV s L,L,a s (Q 2 ) j 



^ [L, a s (Q 2 )] = - ^ f r cusp [ tta (/i 2 )] In |J + 7 H * [a, (/i 2 )] ) , 



(5.8) 



rQ 2 M 



L,L,a s (Q 2 ) 



: Q 2 M 



^ f r cusp [«. (/2 2 )] In ^ - 7 s [a s (/i 2 )] ) , ( 5.0 ) 
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and 



L = In 



1 



L = In ■ 



In {ruQ 2 ) 



(5.10) 



r se< E 

The method of calculating the integral transform in (|5.4|) is developed in Ref . |2jj] . It is 
based on the expansion of the function F s L, L, a s (Q 2 ) into the power series with respect 
to L and using the simple formula 



2tti Jc dg k T(l-g) 



(5.11) 



Using this method, we obtain the following result for the resummed distribution: 

G (r) = [1 + C 3 a s (Q 2 )] as {TQ P TCF expt^^^+^^O,^^ (g ^ 

where a s = a s (Q 2 ) and 



7T 



r[2-g(L,a a 



g{L,a s 



d_ 



F s L, L, a. 



(5.13) 



L=0 



In the expression ( |5.12| ), we have restored the pre-exponential factor 1 + a s (Q 2 ) C3, where 



C3 = Cf 



31 



+ 12C(3) 



+ C A 



353 
~18 



6C(3) 



50 



The corresponding expression for (r) has the form: 

exp [Jh 2 (L, a s ) + F s (L, 0, a s )] 



F (r) = [1 + C 2 a s ] 



where 



C2 — Cf 



r[l-s(L,a.)] 

2vr 2 ' 
-2+ — 



(5.14) 



(5.15) 



(5.16) 



The SCET result ( 5.15| ) presented in Eq.( |5.15 ) coincides with the result of Ref. [p4j] . 

Now one can compare the leading thrust distribution Q5.15 ) and the power suppressed 
one ( pU2|) : 

G(r) a s (tQ 2 ) _ 2 exp {^ g , [L, a s (Q 2 )] - F H2 [L, a s (Q 2 )] } 

~FW) = ~~^ T [1 + asiQ) " C2)] l-9[L, a .{Q*)] ' 

(5.17) 

Taking into account the explicit form of the exponents ( |5 . 9|) , we obtain 

G ( r ) _ r<(Q) f~-\ Mr) ( 5 , 18 ) 



F(t) 



G(°) (r) e" 



where 



(r) = (0) (o) ln (1 _ A) _ ln [x _ g (L) as)] + ^ Cg _ Ca) ^ 

Po 



(5.19) 
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with 



9 (L, a s 



2Tq 



[In (1 — 2A) — In (1 — A)] 



Air 



In 



(5.20) 



Here the pre-exponential factor 1 + a s (Q 2 ) (C3 — C2) has also been exponentiated to the 
NLL level of accuracy. Since the soft factor in Eq.( 4.15| ) is the same as the one in Eq.([4.5|), 
it drops out of the ratio ( |5. 17 ) almost completely. The resummation factor exp [w(r)] is 
shown in Fig. ||. 

Since the function G{t) gives the shape of 
the jets which have distinct angular distribution 
~ sin 2 #3", which is different from that of F(t), 
there is a possibility to measure G(t). Such analy- 
sis was performed by the OPAL collaboration |25| . 
The comparison of theoretical predictions with the 
OPAL data is presented in Fig.|6|, where three curves 
are shown. The thin solid line corresponds to the 
perturbative result (|2.14| ). As one can see from 
Eq.( 3.15|) , dG^/dr tends to constant in the r — > 
limit. The distribution improved by resummation 
( 5.12 ) is drawn by the dashed line. The solid line present the prediction with nonpertur- 
bative effects as discussed below in Sect.||. 

Since the result ( |5.12| ) is valid in the region r <C 1, we match the resummation factors 
with the perturbative result ( 2.14 ) so that all higher order corrections disappear when r 
tends to its maximal value for a three-jet configuration T max = 1/3. The lack of multiplicity 
for the perturbative result ( 2.14 ) explains why the data exceeds the prediction in the region 
r > 1/3, while the poor accuracy of the data in the region r <C 1 does not allow one to 
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Figure 5: Resummation factor as a tunc 
tion of the thrust boundary t > 1 — T. 
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Figure 6: Comparison of the OPAL data for the longitudinal thrust event shape with the theoretical 
predictions dG/dr at r = 1 — T. The bars with short strokes represent systematic errors and those 
with long strokes are statistical errors. The thin solid curve shows O (a s ) perturbative result, the 
dashed curve gives the LL+NLL prediction, while the solid curve is obtained after convoluting with 
the non-perturbative shape function discussed in section [| 
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test the NLL effect specific for G (r) against the background of the common Sudakov 
suppression. 

The factorization scale \j? presented in Eqs.(fPf), ( [4.15 ) formally separates collinear 
intra-jet radiation and soft inter-jet radiation. Since an infrared safe event shape does 
not depend on an explicit definition of jet, the factorization formulae do not depend on 
this scale. This fact is very helpful in deriving the simple representation ( |5.4| ) for G (r) 
where the role of the collinear scale is reduced to the renormalization scale of a s and 
determination of the argument of the hard logarithms in Eq.( |5.17| ). It is worth noticing 
that we use the same r cuS p in the evolution equations ( [5.5| ) for all Hi, although, we test it 



only in the leading order (see Eq.( 3.1$Q ), 



6. Nonperturbative correction 



For the total thrust distribution, the resummation of large logarithms was performed to 
NLL accuracy in Ref . by means of the NLL branching algorithm. The authors of Ref . [ 24 1 
notice that in order to keep the algorithm in the perturbative regime, one has to introduce 
the infrared regulator for the argument of the running coupling constant. Alternatively, 



as it was demonstrated in Refs.[26, 27, p8[ , one may take into account non-perturbative 
(NP) effects by convoluting the resummed perturbative expression with a phenomenological 
shape function u (r): 



NP+PT 



(r) 



rQ/A 



dr'ov 



A 

Q 



U T 



(6.1) 



where A is a phenomenological soft scale characterizing the transition into the NP regime. 
Non-perturbative power corrections generated by this shape function were the subject of 

and theoretical IB4L [3 



intensive experiment al|29|, |30|, |3l], |32 



8|, |36j] studies. Since the 



NP corrections appear mostly due to the radiation of soft partons, NP effects should firstly 
affect the generalized soft factor: 



S T (a,/3;M 2 ) 



£ 

x 



X 



1 n + 1 n 



J) 8(a-n-p x ^j8(p - n+ ■ p Xj ^ . (6.2) 



It implies that NP effects are universal for the distributions containing the same factorized 
soft factor, as confirmed by the analysis performed in Ref. [37]. Since F and G struc- 
ture functions contain the same soft factor, one can estimate the influence of NP effects 
on the ratio G/(F + G) (see Eqs.(4.5) and ( f4.15| )). A very simple but reasonably good 
parametrization for the shape function was found in Ref. 



n K (x) 



r(3/2) 



A = 0.7GeV. 



(6.3) 



Using this parametrization, we find that the NP correction reduces to a simple shift of the 
cross section even in a region r > 1 — (T) = 0.066, so that the relative correction is 



[G/(F + G)] PT , 



NP 



[G/(F + G)] 



PT 



[G/(F + G)] PT 



-0.22 



(6.4) 
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for r = 1 — (T) and increases for smaller r (see Fig. ^ below). 

Since the NP shape function ( |6 . 3| ) is normalized to unity, the NP correction to the 
structure function integrated over total domain, i.e. J drdG/dr, should vanish. However, 
the convolution with uk (x) shifts the point where dG/dr = 0. Thus, if we restrict the 
integration to the true value of r max = 1/2, then the NP correction to G(r max ) would be very 
small. This statement is in agreement with the Monte-Carlo study of hadronization effects 
performed in Ref. ]2| , where some generators found no correction and some generators found 
a small positive correction to the inclusive combination (F — 2G)/(F + 2G). 



7. Forward-backward asymmetry 

In order to suppress QCD corrections to the forward-backward (FB) asymmetry, one can 
reduce final state phase space to the two jet region, thereby suppressing gluon radiation. It 
has been studied in Ref.j|] how the experimental cuts bias the theoretical corrections. The 
event shape can also been used to select the events. Since the phase space of two r-collinear 
partons is of the order r in the r — > limit the corresponding correction decreases with r 
(see Refs.§, |] and Fig.[7|). 

We define the FB asymmetry, which depends on the maximal thrust value (T < 1 — r) , 
as follows 

A / \ _ jo d cos ° T w r ) ~ /-i d cos ° T w ( 6t ' r ) _ ^(0) K ( T ) (71) 
f\dcose T w(6 t ,t) F(t) + G(t)' 

Here 

w(9,r)= f dT da (7.2) 
Ji- T dcost^TdT 

with da (r) defined in Eq.([2.1|), and is the tree-level asymmetry at Q 2 = M% : 

*(0) _ ^ ^9al9vl ^gaq9vq „x 

The effect of QCD radiative corrections can be characterized by the following quantity: 

C{t) = 1-^ = C^- k Ht) + C^(t), (7.4) 
C iF-K) (T) _ F(t)-K(t) 

L ^-F{t) + G{tY (7 ' 5j 



The coefficients C^ - ^) and vanish if the free parton model relations ( 2.11| ), ( |2.12j ) 
hold. The sum F(t) + G(r) gives the integrated thrust cross section. 

We show in Fig. QCD predictions for the correction factor C (r) which relates as 
in Eq.(^J) the electroweak parameter and the observable FB asymmetry of qq jets whose 
thrust value is greater than T > 1 — r. Shown by thin solid lines are the tree- level results 
for C <yF ~ K ' > and C^ G \ The correction factor C^ F ~ K ^ arises from the three parton (qqg) 



-20- 




Figure 7: Correction factors for the forward-backward asymmetry. 



configuration where both q and q are in the same hemisphere along the thrust axis, as 
illustrated in Fig. ||(c). Those events cannot contribute to the P-odd function K (r) and 
the observed asymmetry reduces by the factor C^ F ~ K \ As one can see from Eq. fljOJ ), 
the derivative of the function Fi ree (r) — -ftTtree ( r ) has logarithmic singularity in the two-jet 
limit T — > 1 (see Eq. ( |2,15| )), which prevents us from introducing a local SCET operator 
for the function F — K. Because this region may be studied by parton shower models [38] 
that respect the exact matrix element |3j|] and because the contribution from the relevant 
jet configuration can in principle be removed by double-tag experiments, we concentrate 
our attention on the study of the correction factor in this report. 

In contrast to F — K, the function G (r) originates from the local three-body operator 
3 of Eq. (FTlC| ) in SCET, and we could show in the previous section that the leading 
soft singularities from this operator can be resummed to give the same Sudakov form 
factor as for the total distribution in the two-jet limit (r = 1 — T — > 0). Because this 
common Sudakov factor cancels in the ratio, the thin solid curve for C( G ) in Fig.|7| may 
be regarded as the zeroth-order prediction of perturbative QCD with the leading log (LL) 
thrust distribution. 

It is worth noting here that this interpretation does not hold for the O (a s ) curve for 
(j{F-K) j n pjg ^ because the leading Sudakov factor differs from the qq events even in the 
T — V limit. 

In Fig. [?], our prediction for C^ G ' in the NLL level is shown by the dashed curve 
and thick solid curve is obtained after incorporating the non-perturbative effects by using 
the shape function ( |6.3[ ). Since the NLL level resummation supplemented by the non- 
perturbative correction fl6.3p reproduces the observed thrust distribution rather well (see 
Refs. [p7[ |36[| and Fig. || below) , one can directly compare the FB asymmetry of the observed 
jets in the two-jet region 1 — T < r 1 with the QCD prediction of Eq. ( [7.4[) with the 
correction factor (r) depicted by the thick solid curve in Fig.|7|, provided that the 
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contribution from F — K is reliably excluded by double-tagging and that the r-dependance 
of the observed asymmetry is consistent with the prediction. If one uses all the events with 
T > (T) = 0.934 at Q = mz the QCD correction factor is estimated as 



C (G) 



= 0.0024 ± 0.0002, (7.7) 

=1-(T)=0.066 



in the massless quark limit, when the error is estimated by a quadratic sum of the uncer- 
tainty in a s (Mz) = 0.1184 ± 0.0007 j|l]] and the variation due to the different ansatze for 
the non-perturbative shape function found in Ref. |$7| and Ref. [42|. 

In actual experimental analysis, one needs to correct for effects due to hadronization, 
detector acceptance and performances. They have been accounted for by making use of 
parton shower based event simulator, such as JETSET f38|| , which has been tuned to repro- 
duce all the e + e~ — > hadrons data at collision c. m. energies between 30 and 200 GeV. The 
difficulty of the experimental analysis summarized in Ref. [|| for estimating the electroweak 
asymmetry parameter A^ ' from the observed forward-backward asymmetry of charm and 
bottom quark jets may be traced back to the absence of the parton shower that describes 
the longitudinal structure function in Eq. ( |3.9| ). 

It is worth noting here that the traditional scheme |3t| to match the parton shower 
(that incorporates resummation of LL and NLL emissions) and the matrix elements (that 
give quantum mechanical correlations such as angular distributions) cannot simulate jets 
with correct angular distribution at the accuracy level of precision EW measurements. 
This is because in the leading order, the two-jet like events are matched to the two-parton 
matrix element which in e + e~ —> qq contributes only to the transverse structure functions 
with 1 + cos 2 0t and cos 0t angular distributions. It is only for the three or more jet events 
the tree-level matrix elements contribute to the longitudinal structure function. Although 
the longitudinal contributions are power suppressed in the two-jet region, the accuracy 
required by the EW precision measurements may not allow us to neglect their contribution. 
In this work we have shown that they can be described as parton shower originated from 
the local three-body operator in SCET and that the resulting jet structure is similar but 
different from that of the jets from the qq operator. It is a usual consequence of color 
coherence that two collinear partons are indistinguishable for a large angle radiation and 
should be replaced by a single pseudo-particle in a jet clustering algorithm. The unusual 
point specific for power corrections is that the corresponding three-parton matrix element 
should be supplied by only two A q Sudakov form factors and then combined with two 
jet-like parton showers. The short distance nature of longitudinal events requires modify 
the method to combine parton showers with matrix elements. A successful Monte-Carlo 
simulation of the QCD corrections to the FB asymmetry due to longitudinal structure 
function can be carried out only after such modification is performed. 

Before closing this section let us discuss our predictions for the angular distribution 
of jets with a particular thrust value. The differential jet angular distribution can be 
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expressed as 



do- 3 
d cos 9 T dT ~ 4 



a 



(0) 



F' (1 - T) 



i + cos 2 e T 



+-A^K' (1 - T) cos 6 T + G'(1- T) sin 2 9 T 



(7.8) 



in terms of derivatives of the resummed functions at r = 1 — T. Integration over cos 9t 
gives the thrust distribution: 



g = [F' (1 - D + G" (1 - T)} 
and the FB asymmetry of the jet with a thrust value of T is expressed as 

Jo ldcos ^^dT-/°i dcos ^^ dT 



(7.9) 



A FB (T) 



Icos 6» T dT 



AW. 



J_! d cos flr deodar 
IT' (1 - T) 



F' (1 — T) + G" (1 — T) ' 
The correction factor can also be defined for a given T value 



A FB (T) = A® \l-C (T) 



with 



C(T) 
C(G) ( T ) 



d (F-K) (T) + c-(G) (T) , 
F' (1 - T) - F' (1 - T) 



F'(l 



-T) + G> (1 
G' (1 - T ) 



T) ' 



F' (1 — T) + G' (1 — T) ' 



(7.10) 

(7.11) 

(7.12) 
(7.13) 

(7.14) 



As in the case for the asymmetry of jets with T > 1 — r, in the two-jet region where 
the primary quark and antiquark have momenta in the opposite hemisphere (see Figs. |j](a) 
and (b)), we can safely assume that F (r) = K (r) and all the functions receive common 
non-perturbative corrections. 

The thrust distribution ( [7.9| ) is shown as the dashed curve in Fig. ^ from the resummed 
expressions (5.12) and ( 5.15| ). The inclusion of the non-perturbative soft factor fl6,3| ) shift 
the prediction to the solid curve. The LL+NLL prediction supplemented by the non- 
perturbative soft factor of Eq. ( [T^ ) reproduces the observed thrust distribution in e + e~ 
collisions at all energies between 14 and 206 GeV [ 

gives the mean value (T) = 0.934 and the distribution is peaked at T pea ^ 



37], and in particular at y/s = Mz 

0.9794. Al- 



though our evaluation fl3.16| ) of G' (r) can be regarded as a part of the NNLO correction to 
the thrust distribution, we don't observe significant change in the total thrust distribution. 
In Fig. |, we show our prediction for the QCD correction to the thrust axis angular asym- 
metry, (T), for the jets with a particular thrust value T. Again, the dashed curve 
shows our LL+NLL order prediction and the solid curve is obtained after including the 



- 23 - 




non-perturbative soft factor of Eq. ( p. 3D for both F and G functions. The correction fac- 
tor C (T) decreases by the non-perturbative correction because it reduces the longitudinal 
distribution G' (T) more strongly than F' (T) at T ~ 1. We find 

(T) = 0.0008 at T = T pcak = 0.9794, (7.15) 
(T) = 0.008 at T = (T) = 0.934, (7.16) 

from Fig. ^. The correction factor is small but it can be a significant fraction of the error 
of the 6-jet FB asymmetry which is as small as 1.7% on the Z-boson pole 

Two remarks on our predictions, Eqs.(7.7) and ( 7.16| ), are in order here: The first is on 



the validity of the non-perturbative corrections that lead to our predictions, and the second 
is on the additional correction C^~ K > due to events where both quark and antiquark are 
emitted in the same hemisphere along to the thrust axis opposite the gluon jet direction, 
see Fig.[j](c). 

Our prediction for the asymmetry correction factor C( G > (t) and (T) are ob- 
tained under the assumptions that both the transverse and longitudinal functions F (r) 
and G(t), respectively, receive common non-perturbative corrections via the soft factor 
Eq. (EDI) , whose form has been chosen to reproduce the observed jet thrust distribution in 
e + e~ —> hadrons experiments. We believe that it is an excellent approximation in the two- 
jet (T — > 1) limit where collinear quarks and gluons radiates soft gluons and hadronizes 
coherently. Nevertheless, it is clear that the convolution with the shape function ( |6.3| ) 
does not exhaust all non-perturbative effects. Therefore, a careful study by using a shower 
MC program that incorporates the longitudinal radiation function G (r) is desired for a 
quantitative estimate. 

In this report we consider mainly the correction due to the longitudinal function G (r) 
which arises from the three-jet configurations Figs.|l](a), (b) in the next-to-leading order, 
since the contribution to F — K from the configuration Fig. |l](c) can be removed by re- 
quiring quark and antiquark momenta are in opposite hemispheres in double flavor-tag 
experiments. In practice, however, the double tagging condition leads to to the reduction 
of the number of useful events by one order of magnitude and hence to the loss of accuracy 
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of the measurement. We note here that probably the present PS program matched to the 
tree-level matrix element for three partons |39| can correctly account for those events where 
the quark and the antiquark are in the same thrust hemisphere. It is easy to show that the 
matrix element corresponding to the F — K combination for 1 — T <C 1 is singular in the 
{Pq ~Pg) ~~ ^ limit, where p q , p g are the quark and gluon momenta, respectively. Due to this 
singularity, the logarithmically enhanced contribution comes from the phase-space region 
where the energy of the quark is much smaller than that of the aniquark. For a fixed thrust 
value such hierarchy leads to a large opening angle between q and q. This probably implies 
that one can neglect the interference between the radiation emitted along the quark or the 
antiquark momentum directions, i.e., three-jet configurations give a dominant contribution 
to F — K even for 1-T<1. However, the additional study of accompanied radiation is 
needed. Nevertheless, an extreme care is necessary for a reliable estimate of the error due 
to F — K , since it is the dominant source of the correction to the FB asymmetry as shown 
in Fig.]?]. It mimics the primary quark and antiquark jets in single flavor-tag experiments. 
It is not clear to us if all these points have been appropriately accounted for in the error 
analysis presented in Ref. [g]. 

8. Conclusion 

The corrections to the angular distribution of the thrust axis considered in this paper are 
power suppressed with respect to the event shape variable in the two jet region. We identify 
two mechanisms that radiation affects the angular distribution in the event topology where 
the primary quark and antiquark are radiated into the opposite hemispheres. It is found 
that the short-distance process gives the leading contribution to the longitudinal cross 
section. Using SCET, we propose the factorization formula for the longitudinal cross 
section and perform the large logarithm resummation to the next-to-leading logarithmic 
level of accuracy. The factorization formula allows us to study the leading nonperturbative 
corrections to the longitudinal cross section. A part of the QCD corrections to the forward- 
backward asymmetry is the ratio of the longitudinal and total thrust cross sections. We find 
that the resummation and the nonperturbative effects result in additional suppression of 
this ratio (about 30% for the experimental mean value of the thrust). We observe that the 
short-distance nature of the leading correction yields potential problems with the Monte- 
Carlo simulation of the QCD corrections to the forward-backward asymmetry. We present 
estimates for the QCD corrections to the forward-backward asymmetry in the LL+NLL 
level including non-perturbative corrections. Underestimation of such corrections may be 
relevant for the discrepancy between the weak mixing parameter sin 2 6\y extracted from 
the jet asymmetry data and the others. However, it can be found out only with a help of 
a new improved Monte-Carlo simulation. 
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